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Abstract 

The spin market model [S. Bornholdt, Int.J.Mod.Phys. C 12 (2001) 667] is extended 
into co-evolutionary version, where strategies of interacting and competitive traders 
are represented by local and global couplings between the nodes of dynamic directed 
stochastic network. The co-evolutionary principles are applied in the frame of Bak - 
Sneppen self-organized dynamics [P. Bak, K. Sneppen, Phys. Rev. Letter 71 (1993) 
4083] that includes the processes of selection and extinction actuated by the local 
(node) fitness. The local fitness is related to orientation of spin agent with respect 
to instant magnetization. The stationary regime characterized by a fat tailed dis- 
tribution of the log-price returns with index a ~ 3.6 (out of the Levy range) is 
identified numerically. The non-trivial consequence of the extremal dynamics is the 
partially power-law decay (an effective exponent varies between —0.3 and —0.6) of 
the autocorrelation function of volatility. Broad-scale network topology with node 
degree distribution characterized by the exponent 7 = 1.8 from the range of social 
networks is obtained. 
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The Ising model is the simplest model of the solid state that is an example 
of the N-body problem. From the perspective of economics this example has 
a great importance because it demonstrates that a basic interaction between 
the spins (agents) can bring a non-trivial collective phenomena. The parallels 
between the fluctuations in the economic and magnetic systems affords an ap- 
plication of the spin models to the market statistics. This motivation has been 
outlined by Cont and Buchaud in Ref.jlJ). There is suggested that the equi- 
librium distribution of the "super-spins" can be used to model the stock price 
fluctuations. The "super-spin" model has been extended by the Chowdhury 



* Corresponding author. 

Email address: horvath.denis@gmail.com (Denis Horvath ). 



Physica A 



and Stauffer Ref.f2) by switching on interactions between the lattice spins. 
Further works in this direction Refs. (0; S @) have putted emphasis onto link- 
age of the notion of strategy (as a way of an agent's thinking) with the spin 
degrees of freedom. The combination of these particulars can be made on the 
basis of the famous concept of minority game Ref.(p). The spin variant of mi- 
nority game assumes the competition of the ferromagnetic (short-range) and 
antiferromagnetic (global) couplings that cause an intermittent price dynam- 
ics of the stylised market. The salient points of the model Ref.(|J) are reported 
in below. 

The activity of fundamentalists and interacting traders has been imposed. Each 
interacting trader is characterized by a regular lattice site position i and by the 
corresponding spin variable S^(t) from { — 1, 1}. When S^(i) = 1, the agent 
tends to sell a unit amount of the stock whereas S^'(i) = —1 abbreviates 
the buy order during a given period. The variable S^(i) is updated by an 
asynchronous heat-bath dynamics 



1 with probability P{h {t \i)) , 
-1 with probability 1 - P(/i W (z)) 



where h^'ii) denotes the local field and P(h) is the sigmoid function 

^ = 1+ exp(-2ffl,) < 2 > 



that depends on the inverse fictitious temperature (3. The local field hP^ii) 
reads 

h®(i) = J Yl S®{j)-KS®(i)\m®\. (3) 

j'Gnn(i) 



Here J > is the exchange coupling of nearest neighbors nn(i) and k > is 
the coupling of the local field with instant magnetization 
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where L is the number of the sites at some lattice. When (3 J is bellow the 
critical temperature, the term K,S^(i)\m^\ suppresses the ferromagnetic spin 
phase. For sufficiently large k the frustration between the small (spin) and 
large (lattice) scales yields steady-state bubble-like regime. The determination 
of mutual orientation of spin and minority is only possible with use of integral 
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global information. Mathematically, the negative S^(i)m^ signalizes that 
agent i belongs to the minority. The integral information also affects the price 
of given stock written in the form 



where A is the ratio of the number of fundamentalists to interacting traders. 
The formula can be interpteted as it follows: the predominance of buy or- 
ders implies mS*' > which causes that p^ falls above the fundamental price 
p*. Evidently, the negative rrv® corresponds to under- valued stock. Conse- 
quently, the change m^ t+1 '— rrft' is associated with the logarithmic price return 
ln[p(* +1 )/pW] j n further. 



1 Co-evolutionary dynamics on network of interacting spin agents 

In order to enhance the realism of the model, we have suggested several mod- 
ifications. The primarily inspiration for this originates from the evolutionary 
formulation of minoritv_game that allows agents to adapt strategy according 
to past experience Ref.(|7|). As a result, the co-evolutionary spin market model 
is proposed that imposes the self-organized formation of the stationary dis- 
tributions of strategic variables constituted by the set of inherent fields and 
couplings. 

The static regular lattice geometry represents perhaps one of the most un- 
realistic aspects of aforementioned spin model. To overcome this problem in 
our approach we take advantage of directed network (graph) of labeled nodes 
T = {1,2, . . . ,L}, where node i G T attaches via N ont directed links to its 
neighbors X n (i) 6 T, n = 1, 2, . . . , N ont , i.e. the graph is iV out -regular. The 
links mediate non-symmetric spin-exchange couplings Jjp(i) through which 
the agents may exchange the game-relevant information. Two outgoing links 
X\{i) = 1 + (i)modL, X 2 (i) = 1 + (L + i — 2)modL of node i 6 T are fixed to 
guarantee the network connectedness at any stage of its evolution. The dynam- 
ics of reconnections of the links X n (i), n > 3 (defined below) then develops 
on the background of the quenched bidirectional loop (n = 1,2). In addi- 
tion, the dynamics is constrained to forbid the self-connections (X n (i) = i) 
as well as the multiple connections [X ni (i) = X n2 -t ni {i)]. An important infor- 
mation about the network statistics can be obtained by sampling the distri- 
bution of the node degrees. For our purposes here, we redefine node degree 
fc( in )(j) = Y^=\Hn=i 3j,x n (i) as a number that accounts exceptionally for the 
incoming links of the node j. 

With the above topology in mind, the spin interaction is introduced via mod- 
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Fig. 1. The varying coupling J^\o) shows punctuated equilibrium behaviour, i.e., 
periods of peaceful epochs interrupted by an extremely short fluctuating periods 
of boots and reconnections which are the hallmarks of the self-organized critical- 
ity. We inspected that the probability density function (pdf) of J's indicates the 
formation of the "amorphous" and "glassy" mixture of the couplings that con- 
sists of ferromagnetic (Jn (i) > 0) and (prevailing) antiferromagnetic (J„ (i) < 0) 
contributions. The stationary mean value of J„ (i) equal to —0.518 indicates the 
predominance of antiferromagnetic phase. 

ified local field 

N out 

h®{i) = £ 4 t} (i) S (t) + hg>{i) + «W(<) m« . (6) 

n=l 



We see that the strategic variable /t'''(i) is reminiscent to the previously in- 
troduced global term. The inclusion of the node bias h${i) closely resembles 
a random field model. It should be noted that preliminary simulations have 
confirmed qualitatively that the set of generated distributions is irrespective 
to variations in functional form of the local field. 

As usual, the principles of the co-evolution are simply implementable in the 
frame of the local fitness as a phenotype. Here, the recurrently defined lo- 
cal fitness F®(i) is considered as an integrated history of the agent's spin 
orientations with respect to the minority 

F (t+i) ^ = F (t) ^ _ S (t) ( •) m (t) _ ( 7 ) 

The idea of the slowly (strategic) changing parameters Ref . jij) as a background 
of the spin dynamics is not a new within the family of spin market models. 
In order to optimize the local fitness, the co-evolutionary dynamics including 
the Hebbian and anti-Hebbian learning rule, respectively, has been proposed 
in Ref.(0). The approach proposed here follows similar lines, however, the 
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Fig. 2. The pdf of the log-price returns. In part (a) the index a = 3.6 of N = 10 
links is extracted from fit. The central part of pdf is Gaussian distributed. 

Nearly exponential pdf in part (b) is obtained for N out = 6. We see that pdf 
attained by pruning of links resembles qualitatively Nasdaq composite (c) of 1 min 
data through May - September 2005. 
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The autocorrelation function of the volatility of the log-price returns 
|77iv" ~i — m^'l- The long-time memory with the particular effective decay expo- 
nents is identified numerically. 



Fig. 3. 
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principal distinctions exist that are summarized in below: 

(a) The strategies are subordinated to Bak - Sneppen extremal dynamics 
Refs.JI3;Q|l3); 

(b) In the present formulation we have —S^'m® instead of — [S^ t+1 '(i) — 
S®(i)]m® suggested in Ref.(j^). The effect of the term involving the spin 
difference is that agent benefits from majority /minority jumps, while the 
abidance in minority is paradoxically penalized. 

(c) The dynamical network of spin couplings is permanently rewired which 
mimics the social effects, where emphasis is putted on the links to the fittest 
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Fig. 4. The pdf's of the network node degree k^ n \j) calculated for L = 500, 
_/V out = 10. The comparison with the poissonian distribution obtained for the situa- 
tion when the followers attach to randomly chosen node (no matter how large is the 
fitness). For the attachments to leader and its neighbors imposed the [/c( m )] _7 tail 
is formed. The exponent 7 ~ 1.8 corresponds to the fit within the range < 13, 50 >. 

nodes. Generally, our approach resembles the concept of preferential at- 
tachments that comes from the core of the proposal of the Barabasi and 
Albert Ref.Q). 

Recently, the interest in complex networks has been extended to the search 
for local rules governing to dynamics of social networks. Several principles 
have been exploited for this purpose. As an example let us mention the net- 
work transformation due to node aging and time-dependent attachment prob- 
ability Ref. (jl4l ). The model proposed here is however closer to paradigm of 
inter-agent communication across net Refs. |lil : llfih. 

We have considered slow motion of strategic variables defined by the recursive 
cycles over the following steps: 

(1) The information that leader agent labeled as i max has achieved the highest 
fitness -F(z max ) > F(j), Wj G T is available to all agents. Clearly, the 
idealized and unbounded access to information resources monitoring the 
instant fitness of all agents is assumed. 

(2) The assumption of the model is that the strategic transfer from leader to 
the crowd of L — 1 followers i G T, i 7^ i max is carried out independently of 
the local network structure. The supposition about the collective action 
of the crowd expresses the bounded rationality of the followers. In essence, 
all the followers believe that imitation of strategy that has been applied 
in the aforegoing cycle (i.e. after 1 visit per node) by the most successful 
trader z max would bring a future benefit to them. We postulated that each 
follower i adapts the strategy of leader according to prescription 
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The formula expresses the individual adaption differences. They are re- 
flected by the plasticity parameters r) Jtn (i), r) ka (i), r) R (i) r) PT (i), r] PL {i) 
uniformly distributed over the range (0, r/ max ), where < ?? max < 1. 
(3) What have been left in the previous step without explanation are local 
probabilities of preferential attachments Pi ea d(^), Pt(0- They occur in 
the specific local rules of network dynamics. For given % and randomly 
selected n, n\ > 3 the rewiring follows one of three probabilistic channels 



xJi) 



x n (i) 



is left 

without change 



A-ni (imax) 



with probability 1 — Pi C ad(*) > 

with probability (1 - P T (i)) Pi cad (0 , 
with probability Pi ea d(^) Pr(^) , 



(13) 



where (1 — PT(i))Piead(i) is the probability of the connection to node of 
actual leader. The update X n (i) <— X n , (i max ) favors the social transitiv- 
ity Ref.(|17l). In addition, the transitivity requires bounding Pr(i) < e x ~ 

0(1) > Plead- 

(4) At given t the synchronous update J n (i) <— J^ next/ (i), R(i) <— K^ ncxt/ (i), 

/io(i) <- 4 next) (i), Pi ca d (i) <- Piead Xt) (0 ; PrOO <- Pp (?) is carried out 
for alH G T and n — 1,2, ... , N out . 

(5) The agent i m i n of the lowest fitness P(2 m in) finishes her/his unsuccessful 
life via extremal dynamics. The death-birth rules J^*"*" 1 ^ (*min) — {^ili — 
l)ej, hty%L) = (2r? } - IK, « (t+1) (*2 n ) = (2r? } - l)c« with integer 

^ <+1) («min) = [1 + P^Jjint and Strictly positive P/eld (w) = ^lead, 

vr (*mL) = r 6*'' e T are driven by the uncorrelated uniformly distributed 
numbers {^i*n}n^i\ {^2*n}^=f 5 { r j}}=3 from (0, 1). The last serves for 
"nearly middle" setting p(' +1 ) (i m | n ) = F«(i<l) + r? } (P(4L) -P(^L))- 
The random numbers depend on the free parameters ej, eho, e K , ei ea d, £t 
that bound the available strategic space. The same updates have been 
used for imposition of initial conditions. As one can see from the structure 
of updates, the adaptivity in common with extremal dynamics preserves 
the span of the strategic space. 
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2 Numerical results 



The simulation has been carried out for the spin system thermalized via the 
asynchronous Glauber dynamics with (3 = 1 within the strategic bounds 
eiead = 0.01, ex = 0.5, ej = 8, e K = eh = 4, and with topology L = 500, 
jyout = 10 for t < 2 x 10 7 iterations initialized from random initial settings 
[see step (5) of the algorithm]. As usual, the data corresponding to transient 
regime has been discarded. Hereafter we summarize the main findings of sim- 
ulation. 



The self-organization process of intermittent fluctuations of the log-price re- 
turns has been obtained by the simulation. The punctuated equilibrium of 
selected strategic variable is monitored in Fig 1. 

The relevance of the approach for the economic modeling is demonstrated by 
the fat tail probability density function (pdf) of the log-price returns Ref . (flU) . 
It can be roughly characterized by the exponent a ~ 3.6 (Fig.2). At this 
moment, it is rather illustrative to mention the empirical value a = 3.1 cor- 
responding to 1-day period of TAQ database Ref. ( lib . It has been checked 
numerically that distribution may vary with the bounding of strategic space 
and chosen topology (iV out ). The generalization proposed herein covers not 
only power-law but also nearly exponential or Gaussian-like distributions. 

In Fig. 3 the emergence of the long market memory is demonstrated by calcu- 
lating the autocorrelation of volatility of the log-price returns. The empirical 
slopes —0.34 Ref-lflih and —0.2 Ref. ([2(1) should be mentioned in this context. 
As it turns out from our simulation, the co-evolution represents one of the 
potential mechanisms responsible for the power-law decay. More specifically, 
we have shown the inclusion of extremal dynamics and adaptivity can cure 
the notorious exponential decay problem of the spin market models. Even 
more interestingly, only the simultaneous exploiting of both mentioned dy- 
namical principles can result in nontrivial autocorrelations. Within the spin 
market models, an alternative explanation to our has been proposed only re- 
cently Ref.fjl). In the quoted work, the power-law decay is kept for three-state 
Potts model where zero spin state encodes the inactive (waiting) state of the 
trader. 



In the stationary regime the self-organized reconstruction yields network of 
the broad-scale distributed degrees of nodes (see Fig.4). For the upper part 
of the power-law range we have estimated an exponent 7 ~ 1.8. The re- 
sult is quite encouraging since the exponent falls within the range typical for 
the social networks. Here we report on a comparison the exponents of sev- 
eral cooperative social networks: the value 7 = 1.81 Ref. (|2lh is obtained by 
collecting the e-mal addresses, the exponent 1.2 belongs to the coauthorship 



8 



network Ref. (|22j) whereas 2.1 characterizes the phone-call network Ref. (23). 
Similarly as in Ref. (Eih . the network topology can be characterized by the 
clustering coefficient Cj or its spatial-temporal mean (C). For directed net- 
work we have used slightly modified formula Ci = ej/(A r ° ut (l — iV out )), where 
e; = T,nun^n 3 =i XN ° Ut $x ni (X n2 (i)),x n3 (i) stands for the number of links between 
the neighbors of some node i. In this formula, N out (1 — 7V out ) represents the 
maximum number of links going from X n2 (i) to X n:j (i) and vice versa. It should 
be noted that the conventional factor 2 absents in the definition since the links 
are directed. As usual, it is meaningful to compare the mean clustering coeffi- 
cients of two simulated network regimes. For randomly attached nodes [instead 
of selection declared by Eq.(13)] the simulation results in (C ran d) — 0.027, 
while (C) ~ 0.125 stems from the preference of i max , and thus the enhance- 
ment of clustering is confirmed by the ratio (C) / (C ran d) — 4.6. 



3 Conclusions 



The generalized version of the spin market model has been suggested and 
investigated numerically. In the approach we used the static couplings are ad- 
justed by an instant local fitness that drives the co-evolutionary competitive 
changes. The socially relevant feature of the approach is slowly evolving topol- 
ogy of the spin-spin interaction network as a substrate for fast buy-sell spin 
decisions. The power-law distributions have been uncovered as an attributes 
of the universality caused by the extremal dynamics. Our results suggest hy- 
pothesis that slowly varying strategic variables can cause an emergence of the 
long memory effect that is symptomatic for the real markets. The universality 
manifests itself in a broad-scale node degree distribution, however, the network 
evolution does not necessary imply the power-law decay of volatility. Even the 
static net enough to gain the power-law decay or flat-tailed pdf's. Numerous 
examples of the model applicability can be found outside the field of econo- 
physics. For instance, the coevolutionary model opens the possibility of the 
optimization of the spin lattice and related systems in the manner of Ref. (|25h . 
The finite-size scaling analysis, the small world behaviour and the aspect of 
modularity are only a few examples of issues that remain to be analyzed in 
the complementary studies. 
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